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We present results on lattice QCD at finite temperature (T) and chemical potential {^) . 
We apply the overlap-improving multi-parameter reweighting technique for rif — 4 and 2-1-1 
staggered QCD. We use semi-realistic masses on Lt=4 lattices. The critical endpoint (E) 



of 2+1 flavour QCD on the /^-T plane is located, 
configurations. 



Our results are based on 0(10'' - 10*) 



§1. Introduction 

QCD at finite T and/or fi describes relevant features of particle physics in the 
early universe, in neutron stars and in heavy ion collisions. Extensive experimental 
work has been done with heavy ion collisions at CERN and Brookhaven to explore 
the n-T phase boundary (cf.^^). It is a long-standing question, whether a critical 
point exists on the ^u-T plane, and particularly how to tell its location theoretically.^^ 




Fig. 1. Schematic phase diagram of QCD in the fi — T plane. At /i = and finite T a cross-over is 
expected. The endpoint (E) connects the cross-over region with the first-order region at large 
p. 

Let us start with the /i=0 case first (see Fig. Universality arguments^) and 
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lattice simulations*^) indicate that in a hypothetical QCD with a strange (s) quark 
mass (nis) as small as the up (u) and down (d) quark masses {rriu^d) there would 
be a first order finite T phase transition. The other extreme case (nj=2) with light 
u/d quarks but with an infinitely large rris there would be no phase transition, only 
a crossover. Observables change rapidly during a crossover, but no singularities 
appear. Between the two extremes there is a critical strange mass (rrig) at which 
one has a second order finite T phase transition. Staggered lattice results on Lt=4: 
lattices with two light quarks and nis around the transition T (nj=2+l) indicated^) 
that rrig is about half of the physical m^. Thus, in the real world we probably have 
a crossover. 

At non-vanishing /x, one realizes that arguments based on a variety of models 
(see e.g.^)'^)'^)) predict a first order finite T phase transition at large /i. Combining 
the n = and large fi informations an interesting picture emerges on the /u-T plane. 
For the physical the first order phase transitions at large /j, should be connected 
with the crossover on the /x = axis. This suggests that the phase diagram features 
a critical endpoint E (with chemical potential fj,E and temperature Tg), at which 
the line of first order phase transitions {/i > and T < Te) ends.^^ At this point 
the phase transition is of second order and long wavelength fluctuations appear, 
which results in (see e.g.^-*) consequences, similar to critical opalescence. Passing 
close enough to {iie,Te) one expects simultaneous appearance of signatures which 
exhibit nonmonotonic dependence on the control parameters,^) since one can miss 
the critical point on either sides. 

The location of E is an unambiguous, non-pcrturbativc prediction of QCD. No 
ab initio, lattice analysis based on QCD was done to locate the endpoint. Crude 
models with 7715 = 00 were used (e.g.^)) suggesting that /xg 700 MeV, which 
should be smaller for finite nis- The goal of our work is to propose a new method to 
study lattice QCD at finite and apply it to locate the endpoint. We use full QCD 
with dynamical Uf =2+1 staggered quarks. 

QCD at finite /x can be given on the lattice;^°) however, standard Monte-Carlo 
techniques fail. At Re()u)7^0 the determinant of the Euclidean Dirac operator is 
complex, which spoils any importance sampling method. 

Several suggestions were studied in detail to solve the problem. We list a few of 
them (for a recent review see Ref.^-*^)). 

In the large gauge coupling limit a monomer-dimer algorithm was used.^^) For 
small gauge coupling an attractive approach is the "Glasgow method" in which 
the partition function is expanded in powers of exp(/x/r) by using an ensemble of 
configurations weighted by the /U=0 action. After collecting more than 20 million 
configurations only unphysical results were obtained: a premature onset transition. 
The reason is that the /x=0 ensemble does not overlap sufficiently with the states of 
interest. Another possibility is to separate the absolute value and the phase of the 
fermionic determinant and use the former to generate configurations and the latter 
in observables.^^) The curvature of the /i — T phase diagram at /x = was determined 
by a stochastic calculation of the derivatives of the fermion determinant in.^^) 

At imaginary ^ the measure remains positive and standard Monte-Carlo tech- 
niques apply. The canonical partition function can be obtained by a Fourier trans- 
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form.^^^'^*) In this technique the dominant source of errors is the Fourier transform 
rather than the poor overlap. One can, however, use the fact that the transition hne 
is an analytic function of n, and a polynomial fit for imaginary /i values could be 
analytically continued to real values of fi. The curvature of the phase diagram has 
been determined using this technique for 2,3 and 4 flavour staggered QCD.^^)~^^) 
At temperatures sufficiently above the transition, both real and imaginary n can be 
studied by dimcnsionally reducing QCD.^^^ Hamiltonian formulation may also help 
studying the problem. 



new method 



Glasgow method 



1^ 

Fig. 2. Schematic difference between the present and the Glasgow methods. 



We propose a method to reduce the overlap problem and determine the phase 
diagram in the fi-T plane (for details see^^^). The idea is to produce an ensemble of 

QCD configurations at ;U=0 and at the transition temperature Tc- Then we determine 
the Boltzmann weights^^^ of these configurations at / and at T lowered to 
the transition temperatures at this non-vanishing /x. Since transition configurations 
are reweighted to transition ones a much better overlap can be observed than by 
reweighting pure hadronic configurations to transition ones.^'^) Since the original 
ensemble is collected at fj,=Q we do not expect to be able to decsribe the physics of 
the large /j, region with e.g. exotic colour superconductivity. Fortunately, the typical 
H values at present heavy ion accelerators are smaller than the region we cover. 

Wc apply our technique to 2 + 1 flavour staggered QCD and locate the critical 
point of QCD using semi-physical quark masses. (Multi-dimensional reweighting was 
successful for determining the endpoint of the hot electroweak plasma^^) e.g. on 4D 
lattices.) 
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§2. Overlap improving multi-parameter reweighting 

Let us study a generic system of fermions ip and bosons (j), where the fermion 
Lagrange density is 'ijjM{(p)'ip. Integrating over the Grassmann fields we get: 

Z{a) = j V<PeM-Sbos{oi, </>)] det M(0, a), (2-1) 

where a denotes a set of parameters of the Lagrangian. In the case of staggered 
QCD a consists of /?, ruq and /i. For some choice of the parameters a=ao importance 
sampling can be done (e.g. for Re(/x)=0). Rewriting eq. H2-l() 

Z{a) = j P(/)exp[-S'feos(ao,(A)] det M((/), ao) 

|exp[-S,o,(a,0) + SbosiaoA)]^^^^ • (2-2) 

We treat the curly bracket as an observable (measured on each configuration) and 
the rest as the measure. Changing only one parameter of the ensemble generated 
at ao provides an accurate value for some observables only for high statistics. This 
is ensured by rare fluctuations as the mismatched measure occasionally sampled the 
regions where the integrand is large. This is the overlap problem. Having several 
parameters the set ao can be adjusted to get a better overlap than obtained by 
varying only one parameter. 

The basic idea of the method as applied to dynamical QCD can be summarized 
as follows. We study the system at Re(/i)=0 around its transition point. Using a 
Glasgow-type technique we calculate the determinants for each configuration for a 
set of /i, which, similarly to the Ferrenberg-Swendsen method, can be used for 
reweighting. The average plaquette values can be used to perform an additional 
reweighting in (3. Since transition configurations were reweighted to transition ones 
a much better overlap can be observed than by reweighting pure hadronic configu- 
rations to transition ones as done by the Glasgow-type techniques. The differences 
between the two methods are shown in Figure |21 Moving along the transition line 
was also suggested by Ref.^^^ 

These ideas have been tested using imaginary chemical potential and a perfect 
agreement was found between direct simulations and multiparameter reweighting.-^^) 
Based on these experiences we expect that our method can be applied at Re(/i)7^0. 

§3. The endpoint of = 2 1 QCD 

In QCD with nj staggered quarks one should change the determinants to their 
nj/4 power in our two equations. Importance sampling works also in this case at 
some fi and at Re(/i)=0. Since detM is complex an additional problem arises, one 
should choose among the possible Riemann-sheets of the fractional power in eq. ()2-2() . 
This can be done by using-^^-* the fact that at /i=/iu; the ratio of the determinants is 
1 and it should be a continuous function of /i. 
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Fig. 3. Im(/3o°) as a function of the chemical potential. 



In the following we keep fi real and look for the zeros of Z for complex (5. At 
a first order phase transition the free energy oc logZ(/?) is non-analytic. A phase 
transition appears only in the V— >■ oo limit, but not in a finite V. Nevertheless, Z 
has zeros at finite V, generating the non-analyticity of the free energy, the Lcc-Yang 
zeros. These are at complex parameters (e.g. For a system with first order 
transition these zeros approach the real axis as oo by a 1/F scaling. This V— > oo 
limit generates the non-analyticity of the free energy. For a system with crossover 
Z is analytic, and the zeros do not approach the real axis as oo. 
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Fig. 4. The T-/u diagram. Direct results are given with errorbars. Dotted line shows the crossover, 
solid line the first order transition. The box gives the uncertainties of the endpoint. 

At Tt^O we used I/t=4, ^^=4,6,8 lattices. T=0 runs were done on 10^- 16 lattices. 
"1^,^=0.025 and ms=0.2 were our bare quark masses. At T 7^ we determined the 
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complex valued Lee- Yang zeros, (3q, for different V-s as a function of /i. Their 
CO limit was given by a /3o(^) = Z?^ + C/^ extrapolation. We used 14000, 3600 
and 840 configurations on ^^=4,6 and 8 lattices, respectively. Im(/3^) is shown on 
Figure 131 as a function of /u. For small ^-s the extrapolated Im(/3Q°) is inconsistent 
with a vanishing value, and predicts a crossover. Increasing the value of Im(/3Q°) 
decreases, thus the transition becomes consistent with a first order phase transition 
(overshooting is a finite V effect). Our primary result is l^end — 

0.375(20). 

To set the physical scale we used a weighted average of Rq, nip and ^/a. Note, 
that (including systematics due to finite V) we have {Rq ■ m-,^) = 0.73(6), which is at 
least twice, niu^d is at least four times as large as the physical values. 

Figure m shows the phase diagram in physical units, thus T as a function of fiB, 
the baryonic chemical potential (which is three times larger then the quark chemical 
potential). The endpoint is at Te = 160 ± 3.5 MeV, fiE = 725 ± 35 MeV. At //b=0 
we obtained = 172 ± 3 MeV. 

§4. Conclusion 

We proposed a method - an overlap improving multi-parameter reweighting 
technique - to numerically study non-zero fi and determine the phase diagram in the 
T-n plane. Our method is applicable to any number of Wilson or staggered quarks. 
As a direct test we showed that for lm(^)^0 the predictions of our method are in 
complete agreement with the direct simulations, whereas the Glasgow method suffers 
from the well-known overlap problem. We studied the ii-T phase diagram of QCD 
with dynamical nf=2+\ quarks. Using our method we obtained rE;=160ib3.5 MeV 
and /iE=725±35 MeV for the endpoint. Though is too large to be studied at 
RHIC or LHC, the endpoint would probably move closer to the ^=0 axis when 
the quark masses get reduced. At /i=0 we obtained rc=172ib3 MeV. More work 
is needed to get the final values by extrapolating in the R-algorithm and to the 
thermodynamic, chiral and continuum limits. The details of the presented results 
can be found in.^^^ 

This work was partially supported by grants OTKA-T37615/T34980/T29803/- 
M37071/OMFB1548/OMMU-708 and in part based on the MILC collaboration's 
lattice code: http:/ /physics. india na.edu/^sg/milc. html 

References 

1) P. Braun-Munzinger and J. Stachel, J. Phys. G 28, 1971 (2002). 

2) M. A. Halasz, A. D. Jackson, R. E. Shrock, M. A. Stephanov and J. J. M. Verbaarschot, 
Phys. Rev. D 58, 096007 (1998); J. Berges and K. Rajagopal, Nucl. Phys. B 538, 215 
(1999); M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. Lett. 81, 4816 
(1998); T. M. Schwarz, S. P. Klevansky and G. Papp, Phys. Rev. C 60, 055205 (1999). 

3) R. D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984); F. Wilczek, Int. J. Mod. 
Phys. A 7, 3911 (1992) [Erratum-ibid. A 7, 6951 (1992)]; K. Rajagopal and F. Wilczek, 
Nucl. Phys. B 399, 395 (1993). 

4) For recent reviews see: Z. Fodor, Nucl. Phys. A 715, 319 (2003); E. Laermann and 
O. Philipsen, hep-ph/0303042 S. Muroya, A. Nakamura, C. Nonaka and T. Takaishi, 
Prog. Theor. Phys. 110, 615 (2003); S. D. Katz, hep-lat/0310051 

5) F. R. Brown et at, Phys. Rev. Lett. 65, 2491 (1990); S. Aoki et al. [JLQCD Collaboration], 
Nucl. Phys. Proc. Suppl. 73, 459 (1999). 



Finite T / jjL lattice QCD 



7 



6) A. Barducci, R. Casalbuoni, S. De Curtis, R. Gatto and G. Pettini, Phys. Lett. B 231, 
463 (1989); Phys. Rev. D 41, 1610 (1990); Phys. Rev. D 49, 426 (1994); S. P. Klevansky, 
Rev. Mod. Phys. 64, 649 (1992); M. A. Stephanov, Phys. Rev. Lett. 76, 4472 (1996). 

7) M. G. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B 422, 247 (1998); Nucl. Phys. B 
537, 443 (1999); R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Phys. Rev. Lett. 
81, 53 (1998); K. Rajagopal and F. Wilczek, hep-ph/0011333 

8) S. Borsanyi, A. Patkos, D. Sexty and Z. Szep, Phys. Rev. D 64, 125011 (2001). 

9) M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. D 60, 114028 (1999). 

10) P. Hasenfratz and F. Karsch, Phys. Lett. B 125, 308 (1983); J. B. Kogut et al, Nucl. 
Phys. B 225, 93 (1983). 

11) O. Phihpsen, hep-ph/011005ll 

12) F. Karsch and K. H. Mutter, Nucl. Phys. B 313, 541 (1989). 

13) L M. Barbour et al, Nucl. Phys. Proc. Suppl. 60A, 220 (1998). 

14) D. Toussaint, Nucl. Phys. Proc. Suppl. 17, 248 (1990). 

15) C. R. AUton et al, Phys. Rev. D 66, 074507 (2002). 

16) E. Dagotto et al, Phys. Rev. B41 (1990) 811; 

17) A. Hasenfratz and D. Toussaint, Nucl. Phys. B 371, 539 (1992). 

18) M. G. Alford, A. Kapustin and F. Wilczek, Phys. Rev. D 59, 054502 (1999]^ 

19) P. de Forcrand and O. Philipsen, Nucl. Phys. B 642 290 (2002) arXiv:hep-lat /0205016J. 

20) P. de Forcrand and O. Philipsen, Nucl. Phys. B 673, 170 (2003) arXiv:hep-lat/0307020l. 

21) M. D'Elia and M. P. Lombardo, Phys. Rev. D 67, 014505 (2003) arXiv:hep-lat/0209146 . 

22) A. Hart, M. Laine and O. Philipsen, Nucl. Phys. B 586, 443 (2000); Phys. Lett. B 505, 
141 (2001). 

23) E. B. Gregory, S. H. Guo, H. Kroger and X. Q. Luo, Phys. Rev. D 62, 054508 (2000). 

24) Z. Fodor and S. D. Katz, Phys. Lett. B 534, 87 (2002); JHEP 0203, 014 (2002). 

25) Y. Aoki, F. Csikor, Z. Fodor and A. Ukawa, Phys. Rev. D 60, 013001 (1999); F. Csikor, 
Z. Fodor and J. Heitger, Phys. Rev. Lett. 82, 21 (1999). 

26) A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 63, 1195 (1989); Phys. Rev. Lett. 
61, 2635 (1988). 

27) F. Csikor, G. L Egri, Z. Fodor, S. D. Katz, K. K. Szabo, A. I. Toth, (hep-lat/0104016| 

28) C. N. Yang and T. D. Lee, Phys. Rev. 87, 404 (1952). 



